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This  document  was  prepared  by  the  Electrochemistry  Division,  Directorate  of 
Chemical  Sciences,  Frank  J.  Seller  Research  Laboratory,  United  States  Air 
Force  Academy,  CO.  The  research  was  conducted  under  Project  Work  Unit  number 
2303-F2-10.  Dr.  John  S.  Wilkes  was  the  project  scientist. 
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This  report  has  been  reviewed  by  the  Commander  and  is  releasable  to  the 
National  Technical  Information  Service  (NTIS).  At  NTIS  It  will  be  available 
to  the  general  public.  Including  foreign  nations. 

This  technical  report  has  been  reviewed  and  Is  approved  for  publication. 

C 

CHESTER  J.  Dy«EK,  JR.  Mt  Col,  l£AF 
Director,  Chemical  Sciences 


KENNETH  E.  SIEGENTHALER,  Lt  Col ,  USAF 
Chief  Scientist 

f 

i 

j  Copies  of  this  report  should  not  be  returned  unless  return  Is  required  by 

security  considerations,  contractural  obligations,  or  notice  on  a  specific 
!  document. 

t 

’  Printed  In  the  United  States  of  America.  Qualified  requestors  may  obtain 

;  additional  copies  from  the  Defense  Documentation  Center.  All  others  should 

•  apply  to: 


National  Technical  Information  Service 
6285  Port  Royal  Road 
Springfield  VA  22161 


Unclassified 


security  classification  of  this  page 


1«.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 


2«.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  OECLASSIFICATION/OOWNGRADING  SCHEDULE 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


3.  OISTRI BUT ION/AVAI  LABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
unlimited, 


«.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

FJSRL-TR-86-0004 


5,  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION  |6b.  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 

(If  applicable } 

Frank  J.  Seiler  Research  FJSRL/NC 


JiHIUSjtfSWiTCI 


6c.  ADDRESS  (City,  State  and  ZIP  Coda) 

FJSRL/NC 


7b.  ADDRESS  f City .  State  and  ZIP  Codet 


USAF  Academy 

Colorado  Springs  CO  80840-6528 


8a.  NAME  OF  FUNOING/SPONSORING 
ORGANIZATION 


8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  appticabtei 


AF  Office  of  Scientific  Rsch  AFOSR 


8c.  ADDRESS  (City.  State  and  ZIP  Code I 

AFOSR 
Bldg.  410 


10.  SOURCE  OF  FUNDING  NOS 


PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO. 

NO. 

NO. 

NO. 

°n  6.1 

23Q3 

_ i 

F2 

10 

i2.  personal  AuTHOR<s)Methyl-3-Ethylimid,  Chloride  Rm,  Temp  Molten  Salts 

"11" 


rtftrofTU 


ami 


!lRiTTl 


13a.  type  of  report 


16.  supplementary  notation 


13b.  TIME  COVERED 

from  1  Jun  86  to  30  Jun 


14.  DATE  OF  REPORT  (Yr  .  Mo..  Day /  15.  PAGE  COUNT 

86  Jul  02  40 


COSATl  CODES 


FIELD  GROUP 


18-  SUBJECT  TERMS  (Continue  on  reverte  if  necestary  and  identify  by  block  number) 

density  molten  salts 

viscosity 

conductivity 


19.  ABSTRACT  (Continue  on  reverte  if  necettary  and  identify  by  block  number > 


Density,  specific  conductivity,  equivalent  conductivity,  absolute  viscosity,  kinematic 
viscosity  and  the  Walden  product  for  molten  salts  composed  of  aluminum  chloride  and 
l-methyl-3-ethylimidazolium  chloride  were  plotted  in  three  dimensions,  The  properties 
were  dependent  on  temperature  and  composition.  - 


20  OISTRI8UTION/AVAILABILITY  of  abstract 
UNCLASSIFIED/UNLIMITED  Ijl  SAME  AS  RPT  C  OTIC  USERS  □ 


22a.  NAME  OF  RESPONSIBLE  ING.VIDUAL 

JOHN  S.  WILKES 


DO  FORM  1473,  83  APR 


21  ABSTRACT  SECURITY  CLASSIFICATION 


22b  TE  LEPHONE  NUMBER 
(Include  \  rea  Code ) 


303-472-2655 


EOITION  OF  1  JAN  73  IS  OBSOLETE 


22c  OFFICE  SYMBOL 

FJSRL/NC 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


FJSRL— TR— 86— 0004 


TEMPERATURE  -  COMPOSITION  SURFACES  FOR  PHYSICAL  PROPERTIES 
OF  ALUMINUM  CHLORIDE  -  1-METHYL-3-ETHYLIMIDAZ0LIUM  CHLORIDE 
ROOM  TEMPERATURE  MOLTEN  SALTS 


Lt  Col  Armand  A.  Fannin,  Jr. 
Dr  John  S.  Wi Ikes 


July  1986 


Approved  for  public  release;  distribution  unlimited 


Directorate  of  Chemical  Sciences 
The  Frank  J.  Seller  Research  Laboratory 
Air  Force  Systems  Command 
United  States  Air  Force  Academy 
Colorado  Springs,  Colorado  80840 


LIST  OF  ILLUSTRATIONS 


Figure  1 ,  Density . 4 

Figure  2,  Specific  Conductivity  .  5 

Figure  3,  Equivalent  Conductivity  .  6 

Figure  4,  Kinematic  Viscosity  .  7 

Figure  5,  Log  of  Kinematic  Viscosity . 8 

Figure  6,  Absolute  Viscosity . 9 

Figure  7,  Log  of  Absolute  Viscosity . 10 

Figure  8,  Walden  Product . 11 


> 


INTRODUCTION 

Molten  salts  composed  of  aluminum  chloride  and  l-methyl-3-ethylimidazollum 
chloride  are  Interesting  and  useful  as  electrolytes  in  electrochemical  cells 
and  as  solvents  for  a  variety  of  applications.  We  have  reported  that  the 
physical  properties  (density,  electric  conductivity  and  viscosity)  depend 
strongly  on  temperature  and  melt  composition  (1).  In  our  original  report  we 
plotted  the  physical  properties  as  conventional  two-dimensional  graphs, 
showing  the  dependence  of  temperature  and  compositions  separately.  In  this 
report  we  show  the  properties  plotted  as  three-dimensional  surfaces.  This 
treatment  of  the  data  make  certain  trends  much  easier  to  visualize.  The  plots 
were  made  possible  by  the  acquisition  of  the  graphics  program  DISSPLA”. 
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PLOTTING  METHOD 

DISSPLA*  Is  a  set  of  FORTRAN  calls  that  provides  for  the  plotting  of  data 
or  functional  surfaces  generated  In  user  provided  programs.  Details  about  the 
use  and  content  of  DISSPLA"  are  contained  In  the  software  documentation.  The 
FORTRAN  programs  written  to  plot  figures  1  -  8  are  listed  In  Appendix  A.  The 
essential  feature  of  each  Is  a  function  that  calculates  the  desired  physical 
property  as  a  function  of  temperature  and  composition. 

PHYSICAL  PROPERTY  SURFACES 

The  functions  for  the  temperature  and  composition  dependencies  of  the 
properties  used  were  as  reported  In  reference  1.  In  summary,  the  temperature 
behavior  was  fitted  to  a  V-T-F  model  for  conductivity  and  viscosity  and  to  a 
linear  model  for  density.  The  composition  behavior  for  all  properties  was 
fitted  to  a  complex  structural  model  that  Is  discussed  in  detail  In  reference 
1.  This  model  allows  for  the  marked  change  that  occurs  In  some  properties  at 
equimolar  composition.  The  equations  and  fitted  parameters  are  collected  In 
Appendix  B.  In  all  plots  the  temperature  Is  plotted  In  degrees  Celsius  and 
the  composition  Is  plotted  as  mole  fraction  of  aluminum  chloride. 

Superimposed  on  each  surface  Is  a  polygon  that  encloses  the  area  where  actual 
measurements  were  made. 

Density 

The  density  surface  Is  shown  In  figure  1.  The  density  decreases  linearly 
with  temperature  at  all  compositions,  however  the  strength  of  the  dependence 
varies  with  composition.  The  density  Increases  with  Increasing  A 1  Cl 3  mole 
fraction,  as  expected. 

Electric  Conductivity 

Both  specific  and  equivalent  conductivities  were  plotted.  Specific 
conductivity  Is  shown  In  figure  2  and  equivalent  conductivity  Is  shown  In 
figure  3.  The  cusp  occurring  at  all  temperatures  at  0.5  mole  fraction  Is  much 
sharper  In  specific  conductivity  as  a  consequence  of  the  large  partial  molar 
volume  of  A1C1  . 
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Viscosity 

Kinematic  and  absolute  viscosity  are  presented  In  figures  4-7.  For 
figures  4  and  6  the  units  are  cent1stok.es  and  centlpolse  respectively.  It  may 
be  easily  seen  that  the  melt  viscosity  extends  to  both  higher  and  lower  values 
than  typical  values  for  water.  The  region  of  steep  Increase  in  viscosity 
corresponds  to  the  highly  structured  basic  region  at  low  temperatures.  To 
show  more  detail  In  the  region  of  lower  viscosities,  figures  5  and  7  show 
unitless  representations  where  both  viscosities  are  drawn  on  a  logiQ  scale. 

Walden  product 

For  solutions  where  the  Walden  product  (equivalent  conductivity  times  the 
absolute  viscosity)  Is  a  constant,  changes  In  conductivity  may  be  attributed 
to  only  changes  in  the  viscosity.  While  the  theory  holds  exactly  only  for 
Infinitely  dilute  solutions  In  solvents  of  differing  viscosity,  the  melts  may 
at  least  qualitatively  be  considered  a  solvent  with  viscosity  a  function  of 
the  composition.  Since  concentration  of  the  transporting  species  changes 
smoothly  and  monotonlcal ly  on  each  half  of  the  mole  fraction  diagram,  one 
would  expect  two  monotonlcal ly  varying  Walden  product  surfaces  Intersecting  at 
0.5  mole  fraction.  The  variation  from  such  surfaces  drawn  In  figure  8  gives 
strong  indication  of  some  strong  specific  Interactions  or  structure  In  the 
melts. 
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APPENDIX  A 


The  FORTRAN  programs  for  the  different  three-dimensional  plots  are  largely 
the  same.  The  principal  differences  are  the  function  that  generates  the 
surface,  the  name  of  the  data  file  containing  the  vertices  of  the  polygon 
superimposed  on  the  surface,  the  axes  names  and  ranges  and  the  view  point. 

A  few  comments  about  the  programs  are  Included  here  to  simplify  future  use 
or  modification  for  similar  plots. 

xx(mln)  etc.:  the  minimum  and  maximum  axes  values 

call  ptekal :  Tektronix  4014  terminal  as  output  device.  The  figures 

for  this  report  were  copied  from  a  TEK  4014  screen, 
which  has  the  highest  resolution  available  at  FJSRL  and 
is  much  faster  than  a  digital  plotter, 
call  ptk41 :  Tektronix  4107  as  output  device.  Used  for  ^11  but 

final  copies.  One  or  the  other  of  these  two  calls  Is 
commented  out. 

call  x3name  etc.:  the  axes  labels 

call  vuangl:  establishes  the  viewpoint.  The  arguments  are  with 

respect  to  the  center  of  the  box  containing  the 
surface.  The  first  argument  is  the  angle  (in  the  x-y 
or  temperature-composition  surface)  from  a  line  through 
the  center  of  the  box  containing  the  surface  parallel 
to  the  x  (temperature)  axis.  The  second  argument  is 
the  angle  up  from  the  base  plane  (x,y).  The  third 
argument  is  the  distance  from  the  center, 
call  messag:  Early  versions  had  messages  printed  on  the  box  walls. 

These  were  commented  out  for  the  final  copy, 
call  bshlft:  moves  the  whole  figure  In  x  and  y  directions  In  order 

to  get  complete  figure  on  screen. 


ABSOLUTE  VISCOSITY 
external  avis 
Integer  tic 
logical  ok 

dimension  a<4)  ,b(4)  ,x(2)  ,y(2) ,  laboved ) 
dimension  xx<3) ,yy(3) ,zz(3) 

dimension  xpoly(SO) , ypo1y( 50) ,xcrv( 10) ,ycrv(10) ,zcrv(10) 

data  1above/4HAB0V/ 

mln-1 

max-2 

tic-3 

xx(m1n)-0.0 

xx(max)»100.0 

xx(t1c)-25.0 

yy(m1n)=0.301 

yy(max)-0.70 

yy(t1c)=0. 1 

zz(max)-900. 

zz(m1n)-0.0 

zz(t1c)=300. 

call  1dpoly(xpoly .ypoly.npoly) 
c  call  ptk41 
call  ptekal 
c  call  prpdev(ok) 
c  If (.not.  ok)  stop 

c  call  hp7470 

call  hwrotCAUTO' ) 
call  swlssm 

call  shdchr(90. ,1 ,.002,1) 

call  he  1 gh t ( .325) 

call  physor( .5,0.5) 

call  bshift(0.4,0. ) 

call  nobrdr 

call  area2d<7.5,9.75) 

xl-xmess( 1  Me  1 1  Absolute  VI scosi ty$ * , 100) 

call  messagCMelt  Absolute  Viscosity}'  ,100,4. 0-xl/2.  ,9.25) 

cal  1  volm3d(8. ,8. ,9. ) 

call  mixalf  CINSTRU') 

call  x3name('Temp  (eb.5)o(exbx)C$’ ,100) 

call  y3name('Mole  Fraction  A1C1 ( 1 . 5 ) 3 ( 1 x ) $ ' , 100) 

call  z3name('Abs.  Viscosity  cP$',100) 

call  vuangl(30. ,15. ,40. ) 

cal  1  graf3d(xx(m1 n) ,xx(t1c),xx(max),yy(m1n) , yy <  tic) ,yy(max) , 
1  zz(mln) ,zz(t1c) ,zz(max)) 
call  thkcrv(.05) 
do  25  11-1 .npoly 

cal  1  ntrpol (xpoly, ypoly.npoly , 1 1 , xcrv ,ycrv ,zcrv ,ncrv) 

25  call  curv3d(xcrv,ycrv,zcrv,ncrv,0) 
call  blsur 

call  surfun(av1s,5,xx(t1c)*.2,5,yy(tic)*.2,0) 
call  grfltKO.  ,0.  ,0.  ,0.  ,8.  .0.  ,0.  ,0.  ,9. ) 
call  area2d(8. ,9. ) 
call  cross 

cal  1  graf(yy(m1n) ,yy(t1c) ,yy(max) ,zz(min) , zz<  tic) ,zz(max)) 
call  hei ght< . 525) 


c  x1«xmess< ' AT  Cl  3$ ' ,100) 
c  xpos-4.0-x1/2. 

c  call  messag( * A1 Cl 3$ * ,100,xpos,7.) 

c  call  blrec<xpos-.l ,6.9,xl+.2, .725, .02) 

a(l )»yy(m1n) 
a(2)«yy<max) 
a(3)»yy<max) 
a<4)«yy(m1n) 
b(l )«zz(m1n) 
b(2)«zz(m1n) 
b(3)«zz(max> 
b(4)*zz(max) 
x(l)*yy(m1n) 
x(2)»yy(max) 
y(  1 )»zz(m1n) 
y(2)*zz(m in) 
call  curve(a,b,4,0) 
call  shdpat(19) 
c  shadpat  was  17 

call  shdcrv(x,y,2,0,0,1above) 
call  end3gr(0) 

call  grfltKO.  ,0.,0.,8.,0.,0.,8.,0.,9.) 
call  area2d<8. ,9. ) 
call  cross 

cal  1  graf(xx(m1n) ,xx(tic) ,xx(max) ,zz<m1n) .zz(tlc) ,zz(max)> 
c  next  4  lines  put  msg  on  x-z  wall 
c  xl=xmess< 'deg  C$ ' ,100) 
c  xpos*»4.0-xl /2. 

c  call  messagCdeg  C$' ,100,xpos,7.) 

c  call  blrec<xpos-.l ,6.9,x1+.2, .725, .02) 

a(l )»xx(m1n) 
a(2)»xx(max) 
a(3)=xx(max) 
a(4)=xx(m1n) 
b(l )=zz(m1n) 
b(2)*zz(m1n) 
b(3)=zz(max) 
b(4)-zz(max) 
x( 1 )«xx(m1n) 
x(2)«xx(max) 
y( 1 )«zz(m1n) 
y(2)=zz(m1n) 
call  curve(a,b,4,0) 
call  shdpat(lO) 
c  shadpat  was  18 

call  shdcrv(x,y,2,0,0,1above) 
call  end3gr<0) 

cal  I  grfltKO.  ,0.,0.,8.,0.,0.,0.,8.,0.) 
cal  1  area2d(8. ,8. ) 
call  cross 

cal  1  graf(xx(m1n) ,xx(t1c) ,xx(max) ,yy(m1n) , yy (tic) ,yy(max)) 

a( 1 )-xx(m1n) 

a(2)»xx(max) 

a(3)-xx(max) 

a(4)«xx(m1n) 
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b(l )«yy<m1n) 

b<2)-yy(m1n) 

b(3)«yy(max) 

b(4)-yy(max) 

x(1  )>xx(m1n) 

x(2)-xx(max) 

y(l)«yy(m1n) 

y(2)-yy<m1n) 

call  curve(a,b,4,0) 

call  shdpat<5) 

call  shdcrv(x,y,2,0,0,labove) 

call  end3gr<0> 

call  endpl(O) 

call  donepl 

stop 

end 

function  avls(t.x) 

f-1 . 0-abs  < (2*x-l)/Cl-x)> 

x8»f*f 

If <x. 1 t.0.5>  then 

x2-(1.0-f)*(1.0-f) 
x5=2.0*f*(l  .0-f) 
xl 1-0.0 
X14-0.0 

else 

X2-0.0 

x5-0.0 

x 1 4—  < 1 .0-f )*( 1 .0-f) 
x11=2.0*f*(1.0-f> 

endlf 

a2=252.3 

a5-189.9 

a8=142.2 

all-137.2 

a 14= 132.2 

b2-732 . 5 

b5-0. 

b8-0. 

bl 1-0. 

bl 4-0. 

C2— 4.669 
c5-0. 

C8-0. 
cl  1-0. 
cl  4=0. 

TO»x2*a2+x5*a5+a8*x8+a1 l*xl l+a!4*xl4 
Ck=x2*b2+x5*b5+b8*x8+bl I*x11+bl4*xl4 
aln-x2*c2+x5*c5+c8*x8+cl l*xl 1+c14*x14 
ta-t+273. 15 
ck-732.5 
aln— 4.669 

z»ck/(ta-T0)+0.5*alog(ta)+aln 

elamda-exp(z) 

tb- <  t-60 . 0) * l .0e-4 

z-1 .1279*x2+l .2208*x5+l .2662*x8+l .3263*x1 1+1 .3567*xl4 


dense-z-tb*(6. 1096*x2+6.6170*x5+8.0153*x8+8.7004*xl 1 +9. 61 90* x 14) 
em-146. 64+133. 34*<x/d .0-x)) 

If  (elamda.gt.900. )  elamda-900. 

avl s-(elamda) 

end 

subroutine  ldpo1y(xpoly,ypoly .npoly) 
dimension  xpolyd) , ypolyd) 

openCunlt-S.name-'SdvIs.dat' .status-'OLD' ,err=99) 
do  10  1-1,100 

read<8,  ‘  < 2f  1 2 . 0) '  ,err=99,end=20)  xpolyd )  .ypolyd  ) 
wr1te<5, 1  (2x,2f12.4) 1 )  xpolyd  )  .ypolyd  ) 

10  npoly=1 

20  wr1te(5, ' (2x,14, 1 '  vertex  polygon  from  3dv1s.dat  npoly 

close(8) 
return 

99  wr1te(5,'d'  Error  in  reading  3Dvis.DAT  file*')’) 
stop 
end 

subroutine  ntrpol (xpoly.ypoly , npoly , i i ,xcrv,ycrv,zcrv,ncrv) 
dimension  xpoly<*) .ypolyd) , xcrv(*> .ycrvd) ,zcrv(*) 
nsegs-5 

del ta=l .0/nsegs 
ifdl  .ne.  1)  then 
xl-xpolyd  1-1 ) 
yl-ypolyd  1-1 ) 

else 

xUxpoly(npoly) 
yl=ypoly<npoly) 
end  if 
x2-xpoly( 1 1 ) 
y2=ypoly( i 1 ) 
xx=xl 

del tax=<  x2-xl ) *de 1 ta 

yy=yi 

del tay=(y2-yl )*del ta 
do  30  i-1 ,nsegs+l 
xcrvd  )=xx 
ycrvd  )=yy 
zcrvd  )-avi  s(xx,yy) 
xx-xx+del tax 
30  yy=*yy+del  tay 
ncrv-nsegs+1 
return 
end 


SPECIFIC  CONDUCTIVITY 
external  cond 
Integer  tic 
logical  ok 

dimension  a(4),b(4),x(2),y(2),1above(1) 
dimension  xx(3),yy(3),zz<3) 

dimension  xpo1y(50) ,ypoly(50) ,xcrv< 10) ,ycrv(10> ,zcrv(10) 

data  1above/4HAB0V/ 

mln-1 

max-2 

tic-3 

xx(m1n)-0.0 

xx(max)-120.0 

xx(t1c)-30.0 

yy(m1n)»0.0 

yy(max)-0.70 

yy(t1c)-0. 15 

zz(max)-0.08 

zz(m1n)-0.0 

zz(t1c)-.02 

call  ldpoly(xpoly,ypoly,npoly) 
c  call  ptk41 
call  ptekal 
c  call  prpdev(ok) 
c  If (.not.  ok)  stop 

c  call  hp7470 

call  hwrot('AUTO') 
call  swlssm 

call  shdchr(90.,1,. 002,1) 
call  he1ght( .325) 
call  physor( .5,0.5) 
call  bshlft(0.3,0. ) 
call  nobrdr 
call  area2d(7.5,9.75) 

xl-xmess( 'Melt  Specific  Conduct1v1ty$‘ ,100) 

call  messagCMelt  Specific  Conductivity!’  ,100,4. 0-x1/2.  ,9.25) 

call  volm3d(8. ,8. ,9. ) 

call  mlxalf  CINSTRU') 

call  x3name('Temp  <eh.5)o(exhx)C$' ,100) 

call  y3name('Mole  Fraction  A1C1C1 -5)3<lx)$’ , 100) 

call  z3name('Sp.  Cond.  l/ohm-cm!' ,100) 

call  vuangl (-13. ,62. ,40. ) 

cal  1  graf3d(xx(m1n),xx(t1c) ,xx(max) ,yy(m1n) ,yy(t1c) ,yy(max) , 

1  zz(mln) ,zz(t1c) ,zz(max)) 
call  thkcrv(.05) 
do  25  11-1 .npoly 

call  ntrpoKxpoly, ypo1y.npoly.il  ,xcrv,ycrv,zcrv,ncrv) 

25  call  curv3d(xcrv,ycrv,zcrv,ncrv,0) 
call  blsur 

call  surfun(cond,5,xx(t1c)*.2,5,yy(t1c)*.2,0) 
call  grfltKO.,0.  ,0.,0.,8.,0.,0.,0.,9.) 
call  area2d(8. ,9.) 
call  cross 

cal  1  graf(yy(m1n) ,yy(t1c) ,yy(max) tzz(m1n) ,zz(t1c) ,zz(max)) 
call  he1ght( .525) 
c  x1-xmess( 'A1 Cl  3$' ,100) 
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c  xpos»4.0-x1/2. 

c  call  messag('AlC13$' ,100,xpos,7.) 

c  call  blrec(xpos-. 1 ,6.9,xl+.2, .725, .02) 

a< 1 )»yy(m1n) 
a(2)-yy(max) 
a(3)»yy(max) 
a(4)«yy(mln) 
b(1 )»zz(m1n) 
b(2)»zz(m1n) 
b(3)-zz(max> 
b(4)-zz(max) 
x( 1 )-yy(m1 n) 
x(2)-yy(max) 
y(l )-zz(m1n) 
y(2)-zz(m1n; 
call  curve(a,b,4,0) 
call  shdpat<19) 
c  shadpat  was  17 

call  shdcrv(x,y,2,0,0,1above) 
call  end3gr(0) 

call  grfltKO.  ,8.  ,0.  ,8.  ,8.  ,0.  ,0.  ,8., 9.) 
call  area2d(8. ,9. ) 
call  cross 

cal  1  graf(xx(m1n) ,xx<t1c) ,xx(max) ,zz(m1n) ,zz(t1c) ,zz(max)) 
c  next  4  lines  put  msg  on  x-z  wall 
c  xUxmess(  'deg  C$ ' ,100) 
c  xpos*4.0-xl/2. 

c  call  messagCdeg  C$ '  ,100,xpos,7.) 
c  call  blrec(xpos-. 1 ,6.9, xl+. 2, . 725, .02) 
a< 1 )=xx(m1n) 
a(2)=xx(max) 
a(3)*xx(max) 
a(4)=xx(mln) 
b(l )*zz(m1n) 
b(2)=zz(m1n) 
b<3)=zz(max) 
b(4)*zz(max) 
x( 1 )«xx(min) 
x(2)»xx(max) 
y(  1  )=.zz(m1n) 
y(2)=zz(min) 
call  curve(a,b,4,0) 
call  shdpat(lO) 
c  shadpat  was  18 

call  shdcrv(x,y,2,0,0,1above) 
call  end3gr(0) 

call  grfltKO.  ,0.,0.,8.,0.,0.,0.,8.  ,0.) 
call  area2d(8. ,8. ) 
call  cross 

cal  1  graf (xx(mln) , xx(tlc) ,xx(max) ,yy(m1n) ,yy(t1c) ,yy(max)) 

a( 1 )-xx(min) 

a(2)-xx(max) 

a(3)-xx(max) 

a(4)-xx<m1n) 

b( 1 )-yy(m1n) 
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b(2)-yy(m1n) 

b<3)-yy(max) 

b<4)-yy(max) 

x  ( 1  )-xx(m1n) 

x(2)«xx<max) 

y(l)»yy(m1n) 

y(2)-yy(m1n) 

call  curve<a,b,4,0> 

call  shdpat(5) 

call  shdcrv(x,y,2,0,0, labove) 

call  end3gr(0) 

call  endpl(O) 

call  donepl 

stop 

end 

function  cond(t.x) 

f-1 .0-abs((2*x-1)/(1-x)) 

x8»f*f 

1f(x.lt.0.5)  then 

x2-( 1 .0-f )*( 1 .0-f ) 
x5-2.0*f*(l .0-f) 
x  1 1  *0 . 0 
x 14*0.0 


x2*0.0 

x5*0.0 

xl4-(l .0-f)*<l .0-f) 
xl 1-2 .0*f*< 1 .0-f) 

end  If 
a2*259.9 
a5*189. 1 
a8-123.4 
all-126.5 
al4*121 .9 
b2*493.4 
b5*763.0 
b8*790.9 
bl 1-802.1 
bl 4-813.3 
C2-8.0711 
C5-8.7437 
C8-8.9753 
cl  1-8.9519 
cl  4-8. 9253 

TO-x2*a2+x5*a5+a8*x8+al 1  *  x 1 l+al4*x!4 
Ck-x2*b2+x5*b5+b8*x8+bl  I*xll+bl4*xl4 
aln-x2*c2+x5*c5+c8*x8+cl 1 *xl l+cl4*xl4 
ta-t+273. 15 

z— ck/(ta-T0)-0.5*a1og(ta)+aln 


elamda-exp(z) 
tb-(t-60.0)*l .0e-4 

Z-l .  1279*x2+l .2208*x5+l . 2662*x8+l .  3263*xl 1  +  1 .  3567*xl4 
dense-z-tb*(6. 1096* x2+6. 6170* x 5+8. 0153* x8+8. 7004* xl l+9.6190*xl4) 
em-146. 64+133. 34*(x/<1 .0-x)) 
cond-el amda*dense/em 
end 
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subroutine  ldpo1y(xpoly,ypoly, npoly) 
dimension  xpolyd)  .ypolyd) 

open<un1t-8, name-' 3dC0N.dat' ,status='OLD' ,err»99) 
do  10  1-1,100 

read(8, ' (2fl2.0) '  ,err-99, end-20)  xpolyd )  .ypolyd ) 
wrl te<5, ' <2x,2f 12.4) ’ )  xpoly<  1 ) .ypolyd) 

10  npoly-1 

20  wr1te(5, '  (2x,14, ' '  vertex  polygon  from  3dC0N.dat  ")')  npoly 
c1ose(8) 
return 

99  wr1te<5,'<"  Error  In  reading  3DC0N.DAT  file")') 
stop 
end 

subroutine  ntrpol (xpoly.ypol y, npoly, 1 1 ,xcrv,ycrv,zcrv,ncrv) 
dimension  xpolyd)  ,ypo1y(*)  ,xcrv(*)  ,ycrv(*>  ,zcrv(*) 
nsegs=5 

del ta-1 .0/nsegs 
1 f < 1 1  .ne.  1)  then 
xl=xpoly( 11-1) 
yl=ypolyd  1-1 ) 

el  se 

xl=xpoly(npoly) 
yl=ypoly(npoly) 
end  If 
x2=xpolyd  1 ) 
y2-ypolyd  1 ) 
xx=xl 

deltax-(x2-xl)*delta 

yy=yl 

deltay=(y2-yl)i,delta 
do  30  1=1 ,nsegs+l 
xcrvd  )=xx 
ycrvd  )=yy 
zcrvd  )=cond(xx,yy) 
xx=xx+del tax 
30  yy=yy+deltay 
ncrv=nsegs+l 
return 
end 


DENSITY 

external  dense 
Integer  tic 
logical  ok 

dimension  a(4) ,b(4) ,x(2) ,y(2) , 1  above! 1 ) 
dimension  xx(3) ,yy(3) ,zz(3) 

dimension  xpoly(50) ,ypoly(50) ,xcrv(10) ,ycrv(10) ,zcrv(10) 

data  1  above /4HAB0V/ 

mln-1 

max-2 

tic-3 

xx(m1n)=5.0 
xx(max)=95.0 
xx(t1c)*30.0 
yy(m1n)-0.0 
yy(max)=0.70 
yy(t1c)=0. 15 
zz(max)-1 .50 
zz(m1 n)=l .00 
zz(tlc)-.  1 

call  ldpoly<xpoly,ypoly,npoly) 

C  call  ptk41 
call  ptekal 
c  call  prpdev(ok) 
c  If (.not.  ok)  stop 

c  call  hp7470 

call  hwrotCAUTO' ) 
call  swlssm 

call  shdchr(90. .1 ,.002,1) 

call  he1ght( .325) 

call  physor( .5, .625) 

call  area2d(7. 5,9.75) 

xl-xmess('Melt  Density?' ,100) 

call  messagCMelt  Density?' ,100,4. 0-xl/2.  ,9.25) 

call  volm3d(8. ,8. ,9. ) 

call  mlxalf  CINSTRU') 

call  x3name('Temp  (eh.5)o(exhx)C?' , 100) 

call  y3name('Mole  Fraction  A1CK1 .5)3(1  x)$ '  ,100) 

call  z3name( 'Density  g/mL' ,12) 

call  vuangl (-30. ,30. ,30. ) 

cal  1  graf3d(xx(m1n) ,xx(t1c) ,xx(max) ,yy(m1n) , yy (tic) ,yy(max) , 
1  zz(mln) ,zz(t1c) ,zz(max)) 
call  thkcrv( .05) 
do  25  11-1 .npoly 

cal  1  ntrpol (xpoly.ypoly, npoly, 1 1 .xcrv.ycrv  .zcrv.ncrv) 

25  call  curv3d(xcrv,ycrv,zcrv,ncrv,0) 
call  blsur 

cal  1  surfun(dense,5,xx(t1c)*.2 ,5,yy(t1c)* .2,0) 
call  grfltKO.  ,0.  ,0.  ,0.  ,8.  ,0.  ,0.  ,0.  ,9. ) 
call  area2d(8. ,9. ) 
call  cross 

cal  1  graf(yy(m1 n) ,yy( tic) ,yy(max) ,zz(m1 n) ,zz(t1c) ,zz(max) ) 
call  height! .525) 
xl-xmess( ' A1C1 3?' ,100) 
xpos-4.0-xl/2. 
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call  message  *  A 1  Cl  3$' , 100,xpos,7. ) 

cal  1  blrec(xpos-. 1 ,6.9,xl  +  .2, .725, .02) 

a(l)=yy(m1n) 

a(2)«yy(max) 

a(3)=»yy(max) 

a(4)=yy(ml n) 

b( 1 )=Z2(mi n) 

b(2)=zz(m1 n) 

b(3)=zz(max) 

b(4)=zz(max) 

x( 1 )=yy(m1 n) 

x(2)=yy(max) 

y( 1 )=zz(mi n) 

y(2)=zz(mi n) 

call  curve(a ,b ,4,0) 

call  shdpatC 1 9) 

shadpat  was  17 

call  shdcrv(x,y,2,0,0,1above) 
call  end3gr(0) 

cal  1  grfi ti (0.,8.,0.,8.,8.,0.,0.,8.,9.) 
call  area2d(8. ,9. ) 
call  cross 

cal  1  graf(xx(min) ,xx(t1c) ,xx(max) ,zz(m1n) ,zz(t1c) ,zz(max)) 
next  4  lines  put  msg  on  x-z  wall 
xl=xmess( 'deg  C$' ,100) 
xpos=4.0-xl /2. 

call  message ' deg  C$ ' , lOO.xpos ,7 . ) 

cal  1  blrec(xpos-. 1 ,6.9,xl  +  .2, .725, .02) 

a(l )=xx(min) 

a(2)=xxemax ) 

ae3)=xxemax) 

ae4)=xxemi n) 

b( 1 )=zz(mi n) 

be2)=zz(mi n) 

b(3)=zzemax ) 

be4)=zz(max ) 

xei  )  =  xxemin) 

x  e  2 )  =  x  x  e  ma  x ) 

ye  1 )  =  zzemi n) 

y(2)=zzemi n) 

call  curveea.b.4,0) 

call  shdpatC 1 0) 

shadpat  was  18 

call  shdcrv(x,y,2,0,0, iabove) 
call  end3gre0) 

call  grfiti (0.,0.,0.,8.,0.,0.,0.,8.,0.) 
cal  1  area2de8 . ,8 . ) 
call  cross 

call  graf e  xx(mi n) ,xx(t1c),xxemax) .yyerni n) ,yy(  ti  c)  ,yy(max ) ) 

ae 1 )*xxemi n) 

ae2)=xxemax) 

ae3)=xxemax) 

a(4)=xx(mi n) 

be  1 )=yyemi n) 

b ( 2 ) =yy e  ml n) 


b(3)=yy(max) 

b(4)=yy<max) 

x( 1 )»xx(m1n) 

x(2)»xx(max) 

y(l )»yy(m1n) 

y<2)-yy(m1n) 

call  curve<a,b,4,0) 

call  shdpat(5) 

call  shdcrv(x,y,2,0,0,1above) 

call  end3gr(0) 

call  endpl(O) 

call  donepl 

stop 

end 

function  dense(t.x) 
f-1 -0-abs((2*x-l )/<l-x)) 
tb=*(  t-60 . 0)  *  1  .Oe-4 
x8-f*f 

1f(x. It. 0.5)  then 

x2=(l . 0— f >  *  < 1 .0-f) 
x5=2.0*f*(l .0-f) 
xl 1=0.0 
x 14=0.0 

el  se 

x2=0.0 

x5=0.0 

xl4=(l .0-f )*< 1 .0-f> 
xll-2.0*f*(l .0-f) 

endlf 

Z=1 .1279*x2+l . 2208*  x  5+ 1 .2662*x8+l . 3263*  x 1 1  +  1 . 3567*  x 1 4 
dense=z-tb*(6. 1096*x2+6.6170*x5+8.0153*x8+8.7004*xl 1 +9. 61 90* x 14) 
end 

subroutine  ldpoly(xpo1y,ypo1y,npoly) 
dimension  xpolyd)  ,ypoly(*) 

open(unit=8,name='3dpoly.dat' ,status='0LD' ,err=99) 
do  10  1=1,100 

read(8, ' (2f 12.0) 1 ,err=99,end=20)  xpolyd ) ,ypoly(1 ) 
wr1te(5, '(2x,2fl2.4) 1 )  xpolyd ),ypoly(1 > 

10  npoly=1 

20  wr1te<5, '<2x,14, ’ '  vertex  polygon  from  3dpoly.dat  ")')  npoly 
close(8) 
return 

99  wr1te(5,,(l‘  Error  In  reading  3DP0LY.DAT  file'1)') 
stop 
end 

subroutine  ntrpol (xpoly.ypoly, npoly, 1 1 ,xcrv,ycrv,zcrv,ncrv) 
dimension  xpolyd)  .ypolyd)  .xcrvd)  ,ycrvd)  .zcrvd) 
nsegs=5 

del ta=l .0/nsegs 
1 f < 1 1  .ne.  1)  then 
x!=xpoly( 11-1) 


y1=ypo1y<1  1-1 ) 

else 

x1=xpoly(npo1y) 
yl=ypoly(npo1y) 
end  If 
x2»xpo1y(1 1 ) 
y2=ypoly<11) 
xx=xl 

del tax=<x2-xl )*del ta 

yy,yj 

deltay=(y2-y1 )*delta 
do  30  1=1 ,nsegs+1 
xcrv( 1 )=xx 
ycrv( 1 )=yy 
zcrv( 1 )=dense(xx ,yy) 
xx=xx+de1 tax 
30  yy=yy+de1tay 
ncrv=nsegs+1 
return 
end 
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EQUIVALENT  CONDUCTIVITY 
external  cond 
Integer  tic 
logical  ok 

dimension  a<4),b<4),x(2),y<2),1above(l) 
dimension  xx(3) ,yy<3> ,zz<3) 

dimension  xpo1y<50) ,ypo1y(50) ,xcrv<10) ,ycrv<10) ,zcrv<10) 

data  1above/4HABOV/ 

mln-1 

max-2 

tic-3 

xx(m1n)-0.0 

xx(max)-120.0 

xx(t1c>-30.0 

yy(m1n)=0.0 

yy(max)=0. 70 

yy(t1 c)-0. 1 5 

zz(max)-20. 

zz(m1n)-0.0 

zz(t1c)-5. 

call  ldpoly(xpoly,ypoly,npoly) 
c  call  ptk41 
call  ptekal 
c  call  prpdev(ok) 
c  If (.not.  ok)  stop 
c  call  hp7470 

call  hwrot( 1  AUTO’ ) 
call  swlssm 

call  shdchr(90. , 1 , .002,1 ) 

call  he  1 gh t < .325) 

call  physor( .5,0.5) 

call  bsh1ft(0.3,0. ) 

call  nobrdr 

call  area2d(7.5,9.75) 

xl=xmess( 'Melt  Equlv.  Conducti vi ty$' , 100) 

call  messagCMelt  Equlv.  Conductivity!' ,100,4. 0-xl/2.  ,9.25) 

call  volm3d(8. ,8. ,9. ) 

call  mlxalf  CINSTRU') 

call  x3name( ' Temp  (eh.5)o(exhx)C$‘ , 100) 

call  y3name( 'Mole  Fraction  A1C1 ( 1 . 5)3( 1 x)$ ' , 100) 

call  z3name(’Eq.  Cond.  cm(eh . 5)2(hxex) /ohm-eq$ ' , 100) 

call  vuangl (-1 3. ,62. ,40. ) 

call  graf3d(xx(m1n),xx(t1c) ,xx(max) ,yy(m1n) ,yy(t1c) ,yy(max) , 
1  zz(m1n),zz(t1c),zz(max)) 
call  thkcrv(.05) 
do  25  11-1 ,npoly 

call  ntrpol (xpoly .ypoly.npoly , 1 1 ,xcrv ,ycrv ,zcrv ,ncrv) 

5  call  curv3d(xcrv,ycrv,zcrv,ncrv,0) 
call  blsur 

cal  1  surfun(cond,5,xx(t1c)*.2,5,yy(t1c)*.2,0) 
call  grfltKO.  ,0.  ,0.  ,0.  ,8.  ,0.  ,0.  ,0.  ,9. ) 
call  area2d(8. ,9. ) 
call  cross 

call  graf(yy(m1n) ,yy<t1c) ,yy(max) ,zz(m1n) ,zz(tic)  ,zz(max)) 
call  he  1 gh t < .525) 
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c  xUxmess< 1 A 1  Cl 3$ '  ,100) 
c  xpos-4.0-x1 /2. 

c  call  messag<'AlC13$' ,100,xpos,7.) 

c  call  b1rec(xpos-. 1 ,6.9,x1+.2, .725, .02) 

a(l )-yy<m1n) 
a(2)-yy<max) 
a(3)-yy(max) 
a(4)«yy(mln) 
b( 1 )-zz(m1 n) 
b(2)«zz(m1n) 


b(3)-zz(max) 
b(4)=zz(max) 
x  ( 1  )=*yy(m1  n) 
x(2)=yy(max) 
y( 1 )=zz(mi n) 
y(2)=zz(mi n) 
call  curve(a,b,4,0) 
call  shdpat(19) 
c  shadpat  was  17 

call  shdcrv(x,y,2,0,0,1above) 
call  end3gr(0) 

call  grfltKO.  ,8.  ,0.  ,8.  ,8.  ,0.  ,0.  ,8.  ,9.) 
call  area2d(8. ,9. ) 
call  cross 

cal  1  graf(xx(mln) ,xx<t1c) , xx(max) ,zz(m1n) ,zz(tic) ,zz(max)) 
c  next  4  lines  put  msg  on  x-z  wall 
c  xl=xmess( ' deg  C$' ,100) 
c  xpos=4.0-xl /2. 

c  call  messag( 'deg  C$' , lOO.xpos , 7. ) 
c  call  bl rec( xpos- . 1 ,6 . 9 , xl+ . 2 , . 725, .02) 
a( 1 )=xx(m1 n) 
a(2)=xx(max) 
a( 3)=xx(max) 
a(4)=xx(min) 
b( 1 )=zz(ml n) 
b(2)=zz(mi n) 
b( 3)=zz(max ) 
b(4)=zz(max) 
x(  1  )=xx(fn1n) 
x(2)=xx(max) 
y( 1 )=zz(ml n) 
y(2)-zz(m1n) 
call  curve(a,b,4,0) 
call  shdpatOO) 
c  shadpat  was  18 

call  shdcrv(x,y,2,0,0,1above) 
call  end3gr(0) 

call  grfltKO.  ,0 .  ,0 .  ,  8 .  ,0 .  ,0 .  ,0 .  .  8  .  ,0. ) 
cal  1  area2d(8. ,8. ) 


call  cross 

call  graf(xx(m1n) , xx(tlc) ,  xx(max) ,yy(mln) ,yy(t1c) ,yy(max)) 

a( 1 )-xx(ml n) 

a(2)-xx(max) 

a ( 3 ) - x x ( ma x ) 

a(4)-xx(m1n) 
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b( 1 )-yy(m1n) 

b<2)-yy(m1n) 

b(3)-yy<max) 

b(4)-yy(max) 

x(1 )-xx(m1n) 

x(2)-xx(max) 

y<1  )-yy(m1n) 

y(2)-yy(m1 n) 

call  curve(a,b,4,0) 

call  shdpat(5) 

call  $hdcrv(x,y,2,0,0,1above) 

call  end3gr(0) 

call  endpl(O) 

call  donepl 

stop 

end 

function  cond(t,x) 

f-1 .0-abs((2*x-l )/(l-x)> 

x8-f*f 

if(x.lt.0.5)  then 

x2-(l . 0-f )* ( 1 .0-f) 
x5-2.0*f*(1 .0-f) 
xl 1-0.0 
xl 4*0.0 

else 

x2*0.0 

x5*0.0 

xl 4-( 1 . 0-f ) * ( 1 .0-f) 
x 1 1 —2 . 0*f * ( 1 .0-f) 

endlf 
a2*259.9 
a5- 1 89 . 1 
a8- 123.4 
all-126.5 
a  1 4- 121.9 
b2*493 . 4 
b5— 7  63.0 
b8*790 . 9 
bll-802.1 
bl 4=81 3 . 3 
C2-8.071 1 
C5-8.7437 
C8-8.9753 
cl  1-8.9519 
c 1 4-8 . 9253 

T0-x2*a2+x5*a5+a8*x8+al I*xll+al4*xl4 
Ck-x2*b2+x5*b5+b8*x8+bl l*x1 l+bl4*xl4 
aln«x2*c2+x5*c5+c8*x8+cl l*xl l+c!4*xl4 
ta-t+273 . 1 5 

z— ck/(ta-T0)-0. 5*alog<ta)+aln 

elamda-exp(z) 

tb-<  t-60 . 0) *  1 .Oe-4 

z-1 . 1279*x2+l  .2208*x5+l . 2662*x8+I . 3263*x 1 1 + 1 . 3567#x 14 
dense-z-tb*(6. 1096*x2+6 . 61 70*x5+8 .01 53*x8+8 . 7004* xl 1 +9 . 61 90*xl 4) 
em-146. 64+133. 34*<x/(l . 0- x ) ) 


3 


l  <;»5 


i 


cond-elamda 

end 

subroutine  ldpoly(xpoly,ypo1y, npoly) 
dimension  xpo1y(*) ,ypoly(*) 
open( unit-8,  name-' 3dC0N.dat' .status- 'OLD' ,err-99> 
do  10  1-1,100 

read(8, ' <2fl2.0) ' ,err-99, end-20)  xpoly( 1 ) ,ypoly< 1 ) 
wrl te<5. ' <2x ,2f 12.4) ' )  xpoly( 1 ) ,ypoly( 1 ) 

10  npoly-1 

20  wr1te(5, ‘ < 2 x , 1 4 , ' '  vertex  polygon  from  3dC0N.dat  ")')  npoly 
close(8) 
return 

99  wr 1 te( 5 , ' ( ' '  Error  In  reading  3DC0N.DAT  file'')') 
stop 
end 

subroutine  ntrpol ( xpoly .ypoly .npoly , 1 1 ,xcrv ,ycrv ,zcrv ,ncrv) 
dimension  xpolyC*) ,ypoly(*) ,xcrv(*) ,ycrv(*)  ,zcrv(*) 
nsegs-5 

del t a= 1 .0/nsegs 
1 f < 1 1  .ne.  1)  then 
xl-xpoly(1 1-1 ) 
yl=ypoly( 11-1) 

else 

xl-xpoly(npoly) 
yl-ypoly(npoly) 
end  If 
x2=xpoly(1 1 ) 
y2*ypoly( 1 1 ) 
xx*x  I 

del tax-(x2-xl )*de!ta 
yy-yl 

del tay=(y2-yl ) *  de 1 ta 
do  30  1-1  ,nsegs+ 1 
xcrv( 1 >  =  x  x 
ycrv( 1 )-yy 
zcrv( 1 )-cond( xx ,yy) 
xx-xx+del tax 
30  yy=yy+deltay 
ncrv=nsegs+l 
return 
end 


KINEMATIC  VISCOSITY 
external  avis 
Integer  tic 
logical  ok 

dimension  a(4) ,b(4) ,x(2) ,y<2) , 1above< 1 ) 
dimension  xx(3) ,yy(3) ,zz(3) 

dimension  xpoly(50) ,ypoly(50) ,xcrv(10) ,ycrv<10) ,zcrv(10) 

data  1above/4HABOV/ 

ml  n-1 

max-2 

tl  c-3 

xx(m1 n)-0.0 
xx  (max)- 100.0 
xx(t1c)-25.0 
yy(m1 n)=0. 301 
yy(max)-0. 70 
yy(t1c)-0. 1 
zz(max)-3. 
zz(m1n)-0.0 
zz(t1c)-1 . 

call  ldpoly(xpoly,ypo1y,npoly) 
c  call  ptk41 
call  ptekal 
c  call  prpdev(ok) 
c  If (.not.  ok)  stop 

c  call  hp7470 

call  hwrot( ' AUTO' ) 
call  swfssm 

call  shdchr(90. , 1 , .002 , 1 ) 

call  he  1 gh t ( . 325) 

call  physor( . 5,0. 5) 

call  bsh1ft(0.4,0.) 

call  nobrdr 

call  area2d(7.5.9.75) 

xl»xmess( 'Melt  Kinematic  VI scosl ty$ ' , 100) 

call  messag('Melt  Kinematic  VI scosl ty$ ' , 100,4.0-xl /2 . ,9. 25) 

cal  1  volm3d(8. ,8. .9. ) 

cal  1  ml xalf  ( ' INSTRU' ) 

call  x3name('Temp  (eh.5)o(exhx)C$‘ , 100) 

call  y3name('Mole  Fraction  A 1 C 1 ( 1 . 5)3( 1 x)$' , 100) 

call  z3name('log  Kin.  Viscosity  cS$ ' ,1 00) 

cal  1  vuangl ( 30. , 1 5. ,40. ) 

call  graf 3d( xx(m1 n) , xx(tlc) ,xx(max) ,yy(m1n) ,yy(t1c) ,yy(max) , 
1  zz(m1n),zz(t1c),zz(max)) 
call  thkcrv(.05) 
do  25  11-1 .npoly 

call  ntrpol ( xpoly .ypoly .npoly , 1 1 , xcrv ,ycrv ,zcrv ,ncrv) 

25  call  curv3d( xcrv ,ycrv ,zcrv ,ncrv ,0) 
call  blsur 

call  surfun(av1s,5,xx(t1c)*.2,5,yy(t1c)*.2,0) 
cal  1  grf 1 1 1 (0. ,0. ,0. ,0. ,8. ,0. ,0. ,0. ,9. ) 
call  area2d(8. ,9. ) 
call  cross 

call  graf (yy(m1 n) ,yy( tl c) ,yy(max) ,zz(m1  n) ,zz( 1 1  c) ,zz(max ) ) 
call  he  1 ght ( . 525) 


xUxmess(  ‘ A1C1 3$ ‘  ,100) 
xpos-4.0-xl /2 . 

cal  1  messag( ’ A1 Cl  3$ ’ ,100, xpos, 7. ) 

call  blrec(xpos-. 1 ,6. 9, x1+ .2, . 725, .02) 

a( 1 )-yy(m1n) 

a(2)-yy(max) 

a(3)*yy(max) 

a(4)«yy(m1  n) 

b(l  )«zz(min) 

b(2)-zz(m1 n) 

b(3)»zz(max) 

b(4)*zz(max) 

x( 1 )»yy(m1 n) 

x(2)=yy(max) 

y( 1 )*zz(ml n) 

y(2)=zz(m1 n) 

call  curve(a,b,4,0) 

call  shdpat(19) 

shadpat  was  17 

call  shdcrv(x,y,2,0,0,1above) 
call  end3gr(0) 

call  gr,  itKO.  ,0.  ,0.  ,8.  ,0.  ,0.  .8.  ,0.  ,9.) 
cal  1  area2d(8.  ,9. ) 
call  cross 

call  graf ( xx(mln) ,xx(t1c),xx<max),zz(m1n),zz(t1c) ,zz(max) > 
next  4  lines  put  msg  on  x-z  wall 
xl-xmess( 'deg  C$ ' ,100) 
xpos»4.0-xl /2. 

call  message 'deg  C$ ' , 100, xpos, 7. ) 

call  bl rec( xpos- . 1 ,6.9, xl+. 2, .725, .02) 

a( 1 )=xx(m1 n) 

a  ( 2 )  =  x  x  ( ma  x ) 

a(3)=xx(max ) 

a(4)=xx(m1n) 

b( 1 )=zz(mi n) 

b(2)=zz(min) 

b( 3)=zz(max ) 

b(4)=zz(max) 

x( 1 )=xx(m1 n) 

x (2 ) =xx(max ) 

y(  1 )*zz(min) 

y(2 )  =  zz(mi n ) 

call  curve(a ,b,4,0) 

cal  1  shdpat( 10) 

shadpat  was  18 

call  shdcrv(x,y,2,0,0,1above) 
call  end3gr(0) 

call  grfltKO.  ,0.  ,0.  .8.  .0.  ,0.  ,0.  ,8.  ,0. ) 
call  area2d(8.  ,8. ) 
call  cross 

call  graf(xx(min),xx(t1c),xx(max),yy(m1n),yy(t1c),yy(max)) 

ad)-xx(mln) 

a  ( 2 )  -  x  x  ( ma  x ) 

a  ( 3 )  =•  x  x  ( ma  x ) 

a(4)-xx(m1n) 


b( 1 )-yy(m1 n) 

b(2)-yy(m1n) 

b(3)-yy(max) 

b(4)»yy(max) 

x( 1 )-xx(m1n) 

x(2)-xx(max) 

y( 1 )-yy(ml n) 

y(2)-yy(m1 n) 

call  curve(a,b,4,0) 

call  shdpat(5) 

call  shdcrv(x,yl2,0,0,labove) 

call  end3gr(0) 

call  endpl(O) 

call  donepl 

stop 

end 

function  avls(t.x) 

f-1 .0-abs((2*x-l  )/< l-x> ) 

x8=f*f 

If (x . 1 t .0. 5)  then 

x2=C 1 .0-f)*(l .0-f) 
x5-2.0*f*(l .0-f) 
xll-0.0 
x 14-0.0 

else 

x2=0.0 
x 5=0.0 

xl  4»(  1 .0-f )*( 1 .0-f) 
xl 1 -2 . 0*f *( 1 .0-f) 

endl  f 
a2-252 . 3 
a5-189. 9 
a8-l 42 . 2 
all-137.2 
a  1 4=1 32 . 2 
b2»732 . 5 
b5-0 . 
b8-0 . 
bl 1-0. 
bl 4-0 . 

C2— 4.669 
C5=0. 
c8-0 . 

Cll-O. 
cl  4-0. 

T0-x2*a2+x5*a5+a8*x8+a1 l*xl l+al4*x!4 
Ck.-x2*b2  +  x5*b5+b8*x8+bl  l*x1 l+bl4*xl4 
a1n-x2*c2+x5*c5+c8*x8+cl  l*xl l+cl4*xl4 
ta-t+273. 1  5 
ck-732.5 
a  In— 4.669 

z«ck/( ta-T0)+0. 5*a1og( ta)+a1n 

elamda-exp(z) 

tb-<  t-60 . 0) *  1 .Oe-4 

de-1 . 1 279*x2  +  l . 2208*x5+1 . 2662* x 8+1 . 3263*  x  1  1  +  1  . 3567  *  x  1 4 
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dense-de-tb*(6. 1096* x2+6. 6170* x5+8. 0153* x8+8. 7004* xl l+9.6190*xl4) 

em-146.64+133.34*(x/(l ,0-x)) 

e 1 amda-e 1 amda/dense 

If  (elamda.gt. 1000. )  elamda-1000. 

avls-aloglO(elamda) 

end 

subroutine  ldpoly(xpoly,ypoly, npoly) 
dimension  xpoly<*) ,ypoly(*) 

open ( unit-8,  name- '3dv1s.dat' .status* 'OLD* ,err=99) 
do  10  1-1,100 

read(8, 1 <2f 12.0) 1 ,err-99, end-20)  xpolyd  > ,ypoly( 1 ) 
wrl te(5 , ' (2x ,2f 12 .4) 1 )  xpolyd ) ,ypo1y( 1 ) 

10  npoly=1 

20  wr1te(5, ' (2x,14, 1 '  vertex  polygon  from  3dv1s.dat  npoly 

close(8) 
return 

99  wr i te ( 5 , ' ( * '  Error  In  reading  3Dvis.DAT  file'')') 
stop 
end 

subroutine  ntrpol (xpoly .ypoly, npoly , 1 1 ,xcrv,ycrv,zcrv,ncrv) 
dimension  xpolyd)  ,ypo’  v(*)  ,xcrv(*)  ,ycrv(*)  ,zcrv(*) 
nsegs=5 

delta-1 .0/nsegs 
1 f ( 1 i  .ne.  1)  then 
xl=xpoly( 11-1) 
yl=ypoly( 11-1) 

else 

xl-xpoly(npoly) 
yl-ypoly(npoly) 
end  if 
x2=xpolyd 1 ) 
y2=ypoly( i i ) 
x  x  =  x  1 

deltax  =  (x2-xl  )*delta 
yy=yl 

del tay=(y2-yl )*del ta 
do  30  1=1 , nsegs+1 
xcrv( 1 )  =  x  x 
ycrv( 1 )=yy 
zcrv( 1 )=avi s( xx ,yy) 
xx  =  xx+de 1  tax 
30  yy=yy+deltay 
ncrv=nsegs-t-1 
return 
end 


WALDEN  PRODUCT 
external  wald 
Integer  tic 
logical  ok 

dimension  a(4) ,b(4) ,x<2) ,y(2) , laboved ) 
dimension  xx(3) ,yy(3) ,zz(3> 

dimension  xpo1y<50> ,ypoly(50) ,xcrv(10) ,ycrv(10) ,zcrv(10> 

data  1  above /4HABOV/ 

mln-1 

max-2 

tic-3 

xx(m1 n)=0.0 
xx(max)=«100.0 
xx(t1c)*25.0 
yy(m1n)=0. 301 
yy(max)=0. 70 
yy(t1c)=0. 1 
zz(max)=l 6000. 
zz(m1n)=8000.0 
zz(t1c)=2000. 

c  call  ldpoly(xpoly,ypoly,npoly) 

c  call  ptk41 

call  ptekal 
c  call  prpdev(ok) 
c  If (.not.  ok)  stop 
c  call  hp7470 

call  hwrot( 'AUTO* ) 
call  swlssm 

call  shdchr(90. ,1 ,.002,1) 

call  helghtC .325) 

call  physor( .5,0.5) 

call  bsh1ft(0. 3,0. ) 

call  nobrdr 

call  area2d(7.5,9.75) 

xl-xmess( 'Mel t  Walden  Products* ,100) 

call  messag('Melt  Walden  Products* ,100,4. 0-xl/2. ,9.25) 

call  volm3d(8. ,8. ,9.) 

call  mlxalf  CINSTRU') 

call  x3name('Temp  (eh . 5)o(exhx)C$ ' ,100) 

call  y3name('Mole  Fraction  A1CK1 .5)3(lx)$*  ,100) 

call  z3name( ' (m6)L(mx)  *  (m7)c(mx)$* , 100) 

call  vuangl (-1 1 . ,40. ,40. ) 

cal  1  graf3d(xx(m1n) ,xx(t1c) ,xx(max) ,yy(m1n) , yy( tic) ,yy(max) , 
1  zz(mln) ,zz(t1c) ,zz(max)) 
call  thkcrv(.05) 
do  25  11—1 .npoly 

c  call  ntrpoKxpoly, ypoly, npoly, 11  ,xcrv,ycrv,zcrv,ncrv) 

25  call  curv3d(xcrv,ycrv,zcrv,ncrv,0) 
call  blsur 

call  surfun(wald,5,xx(t1c)*.2,5,yy(t1c)*.2,0) 
call  grfltKO.  ,0.  ,0.  ,0.  ,8.  ,0.  ,0.  ,0.  ,9. ) 
call  area2d(8. ,9. ) 
call  cross 

cal  1  graf(yy(m1n) , yy( tic) ,yy(max) ,zz(m1n) ,zz(t1c) .zz(max)) 
call  he  1 gh t ( . 525) 


c  x1=xmess( ' A1 Cl  3$ ' ,100) 
c  xpos*4.0-x1/2. 

c  call  message ‘ A 1  Cl 3$ ' ,100,xpos,7. ) 
c  call  b1rec(xpos-. 1 ,6.9,x1+.2, .725, .02) 
a(1)-yy<min) 
a<2)*yy(max) 
a(3)*yy(max) 
a(4)»yy(m1n) 
b( 1 )»zz(ml n) 
b(2)*zz(m1n) 
b(3)»zz(max) 
b(4)=zz(max) 
x  C 1  )=yy(min) 
x(2)=yy(max) 
y( 1 )=zz(mi n) 
y(2)=zz(m1n) 
call  curve(a,b,4,0) 
call  shdpat(19) 
c  shadpat  was  17 

call  shdcrv(x,y,2,0,0,1above) 
call  end3gr(0) 

call  grfltKO.  ,8.  ,0.  ,8.  .8.  ,0.  ,0.  ,8.  ,9. ) 
call  area2d(8. ,9. ) 
call  cross 

cal  1  graf(xx(m1n) ,xx(t1c) ,xx(max) ,zz(min) , zz(tlc) ,zz(max)) 
c  next  4  lines  put  msg  on  x-z  wall 
c  xl=xmess( 1  deg  C$' ,100) 
c  xpos=*4.0-xl/2. 

c  call  messag< 'deg  C$' ,100,xpos,7. ) 

c  call  blrec(xpos-. 1 ,6.9,xl+.2, .725, .02) 

a(l )=xx(mln) 
a(2)=xx(max) 
a(3)=xx(max) 
a(4)*xx(mi n) 
b( 1 )=zz(mi n) 
b(2)=zz(m1 n) 
b(3)=»zz(max) 
b(4)=zz(max) 
x( 1 )axx(mi n) 
x(2)*xx(max) 
y( 1 )=zz(mi n) 
y(2)*zz( min) 
call  curve(a,b,4,0) 
call  shdpat(lO) 
c  shadpat  was  18 

call  shdcrvex.y^.O.O.labove) 
call  end3gr(0) 

call  grfltKO.  ,0.  ,0.  ,8.  ,0.  ,0.  ,0.  ,8  .  ,0. ) 
cal  1  area2d(8. ,8. ) 
call  cross 

cal  1  graf (xx(mln) ,xx(t1c),xx(max),yy(m1n),yy(t1c) ,yy(max)) 

a( 1 )*xx(m1 n) 

a(2)=ixx(max) 

a(3)«xx(max) 

a(4)*xx(m1 n) 


b<l )=yy(m1n) 

b(2)«yy<m1n) 

b(3)=yy(max) 

b(4)=yy(max) 

x(1>«xx(m1n) 

x<2)=xx(max) 

y<l )=yy(m1n) 

y(2)=yy(m1n) 

call  curve<a,b,4,0) 

call  shdpat(5) 

call  shdcrv(x,y,2,0,0, labove) 

call  end3gr(0) 

call  endpl(O) 

call  donepl 

stop 

end 

function  wald(t,x) 

f=l .0-abs( (2*x-1 ) /( 1-x) ) 

x8=f‘f 

1fCx.lt. 0.5)  then 

x2=( 1 .0-f )*( 1 .0-f ) 
x5=2.0*f*(l .0-f) 
xl 1=0.0 
xl4=0.0 

else 

x2=0.0 

x5=0.0 

x!4=( 1 .0-f )*< 1 .0-f > 
xl l=2.0*f*(l .0-f) 

endlf 
a2=259.9 
a5=189. 1 
a8=123.4 
all-126.5 
a  1 4= 121.9 
b2=493.4 
b5=763.0 
58=790.9 
bl 1=802.1 
b 1 4=8 13.3 
c2=8.071 1 
c5=8 . 7437 
c8=8. 9753 
cl  1=8.9519 
c 1 4=8 .9253 

T0=x2*a2+x5*a5+a8*x8+al 1 *x1  Hal4*xl4 
Ck=x2*b2+x5*b5+b8*x8+bl I*xl l+b!4*xl4 
aln=x2*c2+x5*c5+c8*x8+cl 1 *xl l+c!4*x14 
ta=t+273. 15 

z=-ck/(ta-T0)-0.5*a1og(ta)+aln 

elamda*exp(z) 

tb=<  t-60 . 0) *  1 .0e-4 

z«l . 1279*x2+l .2208*x5+l .2662*x8+1 . 3263*xl 1+1 .3567*xl4 
dense=z-tb*(6. 1096*x2+6. 61 70*x5+8.01 53*x8+8. 7004*xl 1+9.61 90*xl4) 
em=146. 64+133. 34*<x/(1.0-x>> 


el  amda=elamda*em/dense 

av2=252.3 

av5*189. 9 

av8=142.2 

avl 1=137.2 

av  1 4=.  1 32 . 2 

etak=732.5 

etaaln— 4.669 

etaT0=x2*av2+x5*av5+av8*x8+avl 1 *xl 1+avl4*xl4 
avi s=exp(etak/(ta-T0)+0.5*alog(ta)+etaaln) 
wald=m1n(elamda*av1s, 20000. ) 
end 

subroutine  1dpoly<xpo1y,ypoly .npoly) 
dimension  xpoly(*) ,ypoly(*) 

open(uni t=8,name=' 3dC0N.dat 1 , status* 'OLD' ,err=99> 
do  10  1=1,100 

read(8, 1 (2fl2.0) ' ,err=99,end=20)  xpoly(l) ,ypoly(i ) 
wri te(5, ' (2x ,2f  1 2 .4) ' )  xpol y( i ) ,ypoly( i ) 

10  npoly=i 

20  write(5, ' (2x,i4, 1 '  vertex  polygon  from  3dC0N.dat  ,l)')  npoly 
close(8) 
return 

99  write(5,'(''  Error  in  reading  3DC0N.DAT  file11)') 
stop 
end 

subroutine  ntrpol ( xpol y.ypoly, npoly, 1 1 ,xcrv,ycrv,zcrv,ncrv) 
dimension  xpoly(*) ,ypoly(*) , xcrv(*> ,ycrv(*> ,zcrv(*) 
nsegs=5 

del ta=l .0/nsegs 
if  ( i i  . ne .  1 )  then 
xl=xpoly( i 1-1 ) 
yl=ypoly(i 1-1 ) 

el  se 

xUxpoly(npoly) 
yl=ypoly(npoly) 
end  if 
x2=xpoly( i i ) 
y2=ypoly( 1 1 ) 

X  X  =  X  1 

del tax=( x2-xl ) *  de 1 ta 

yy=yi 

deltay=(y2-yl)*delta 
do  30  i=l ,nsegs+l 
xcrv( i )=xx 
ycrv( 1 )=yy 
zcrv( i )=wald(xx ,yy) 
xx=xx+del tax 
30  yy=yy+deltay 
ncrv=nsegs+l 
return 
end 


APPENDIX  6 


THREE-DIMENSIONAL  EQUATIONS  FOR  MeEtlmCl-AlCl 3  BINARY  MELTS 

DENSITY:  g/cma  [Valid  for  0.25  <  N  <  0.6666  and  10  °C  <  t  <  90  °C] 
p  =  Xa+Xa+Xa+Xa  +  X  a  + 

-2~2  -5~5  -8-8  -11-11  -14-14 

(Xb  +  X  b  +Xb  +X  b  +  X  b  )  (t  -  60)  (1) 

-2-2  -5-5  -8-8  -11-11  -14-14  ~ 

where  the  a^ 1  s  and  fc^'s  represent  the  parameters  In  the  equation  below  for  the 
Individual  "complexes"  with  2,  5,  8,  11,  and  14  chlorides. 

P  =  a  +  b  (t  -  60)  (2) 

The  X . ' s  may  be  calculated  most  conveniently  as  follows:  Let  f  be  the 
anion  fraction  of  A 1 C 1 ^  ~.  Then  f  =  1  -  |(2N  -  1)/(1  -  N) | .  For  0  <  N  <  0.5,  X?  * 
(1  -  f)2,  X  =  2f(l  -  f),  and  X  and  X  =  0.  For  0.5  <  N  <  0.667,  X  and  X  *  0, 
X  =  2f(l  -  f),  and  X  =  (1  -  f)2.  At  all  compositions  X  =  f2. 

-11  -  -  -14  -  -a  - 

TABLE  I,  Least  Squares  Fitted  Parameters  for  Equation  1  Densities. i 


1  *1 


-bj  X  104 


1.1279 
1 .2208 
1.2662 
1.3263 
1.3567 


6.1096 

6.6170 

8.0153 

8.7004 

9.1690 


-  The  parameters  yield  a  »  0.0004  in  p,  where  a  =  [E(pcalc  ~  Pexp>*/ 

1  /2 

(number  of  points  -  number  of  fitted  parameters)] 


37 


CONDUCTIVITY:  cnf/eq/L  [Valid  for  0.30  <  N  <  0.6666  and  20  °C  <  t  <  100  °C] 
VISCOSITY:  cP  [Valid  for  0.31  <  N  <  0.6666  and  10  °C  <  t  <  90  8C] 


and 


where 


ln  A  - 

-kA  < 

(T  -  T  ) 

-  1/2  In  T  +  ln  Aa 

(3) 

-A 

-o 

-A 

1  n  n  = 

k  / 

(T 

-  T  ) 

+  1/2  In  T  +  ln  A 

(4) 

-n 

-o 

~n 

T 

-  X  a 

+ 

X  a  + 

Xa+Xa  +Xa 

(5) 

-o 

-2-2 

-5-5 

-8-8  -11-11  -14-14 

k 

=  X  b 

+ 

X  b  + 

Xb+Xb  +  X  b 

(6) 

-2-2 

-5-5 

-8-8  -11-11  -14-14 

ln  A 

=  X  c 

+ 

X  c  + 

Xc+Xc  +Xc 

(7) 

-2-2  -5-5  -8-8  -11-11  -14-14 

where  the  a’s,  b's,  and  c‘s  represent  the  T  ,  k  and  In  A  respectively. 


-o-  -A  -A’ 

for  each  of  the  Individual  "complexes"  with  2,  5,  8,  11,  and  14  chlorides 


The  X^'s  are  calculated  as  given  for  densities. 

TABLE  II.  Least  Squares  Fitted  Parameters  for  Equations  3,  4  and  5-7  Equivalent 
Conductl vl ties  and  Absolute  Viscosities  of  1 -Me th y 1-3-ethyl imldazol lum  Chloride 
Aluminum  Chloride  Binary  Mixtures 


Equations  5-7  "complex"  species  parameter 


Eq  3  or  4 


Parameters 

2 

5 

8 

11 

14 

Eq  3  and 

5-7  Parameters  for 

Equivalent  Conductivity- 

Io'K 

259.5 

189.1 

123.4 

126.5 

121.9 

!Sa’K 

493.4 

763.0 

790.9 

802.1 

813.3 

1  n  A. 

8.0711 

8.7437 

8.9753 

8.9519 

8.9253 

-a 

Eq  4  and  5-7  Parameters  for  Absolute  Viscosity 

b 

Io'K 

252.3 

189.9 

142.2 

137.2 

132.2 

k  ,K 
— n 

732.5- 

In  A 
~n 

-4.669- 

i  The  parameters  yield  a  -  0.02  In  In  A,  where  a  =*  [Kin  Acaic  -  In  Aexp)2/ 
(number  of  points  -  number  of  fitted  parameters)]1^, 
k  The  parameters  yield  o  =  0.03  in  In  n,  where  a  =*  [Kin  ncalc  ~  ln  nexp)z/ 
(number  of  points  -  number  of  fitted  parameters)]1^. 

£  These  values  are  not  parameters  but  averaged  values  of  kn  and  ln  An  to  be 
used  for  absolute  viscosity  calculations  at  all  compositions. 


